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Abstract 

An application of the particular type of nonlinear operator algebras to spectral problems 
is outlined. These algebras are associated with a set of one-dimensional self-similar potentials, 
arising due to the g-periodic closure fj + pf(x) = qfj(qx), kj + N = q 2 kj of a chain of coupled 
Riccati equations (dressing chain). Such closure describes g-deformation of the finite-gap 
and related potentials. The N = 1 case corresponds to the g-oscillator spectrum generating 
algebra. At iV = 2 one gets a g-conformal quantum mechanics, and N = 3 set of equations 
describes a deformation of the Painleve-IV transcendent. 
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Quantum algebras, or g-deformations of Lie algebras, attracted much attention during the last 
years. They have appeared in conformal field theories, spin-chain models, in the construction 
of link invariants, and so on. Another type of nonlinear algebras, known as W-algebras, also 
formed the subject of recent intensive investigations. It is quite natural to try to find applications 
of such objects within the context of Sturm-Liouville type spectral problems. This is not so 
straightforward and most of the attempts done in this direction cary phenomenological character 
- the deformations of spectra are introduced in an ad hoc manner, without proper definition of 
arising operators. As we shall show below, the problem can be treated in a rigorous fashion on the 
basis of standard concepts of continuous space physics. Moreover, nonlinear algebras will appear 
very naturally as an inevitable consequence of the analysis of exactly solvable potentials. 

The consideration will be limited to the simplest one-dimensional Schrodinger equation 

- ijj'{x) + U{x)i){x) = Xip(x), (1) 

describing a particle moving in some potential U(x). (The prime in ([I]) and below always denotes 
derivative with respect to real coordinate x, x G R.) The quantum mechanical spectral problems, 
associated to ([!]), traditionally provide good place for probing new group-theoretical ideas. This 
is inspired by the general belief that any kind of regularity in spectral data is generated by 
some symmetry algebra of a Hamiltonian. The qualitative understanding of peculiarities of a 
given system is reached when an operator algebra governing the map of physical states onto 
each other is explicitly constructed. The most advanced approach to building of such symmetry 
transformations, which we are going to exhibit, is connected with the technique developed by 
Darboux long time ago IJ. Within quantum mechanics, it is known as factorization method 
[J2J. In the theory of integrable nonlinear evolution equations, it was generalized and named 
as the dressing method ||. Later, Darboux transformations were identified as supersymmetry 
transformations, mixing bosonic and fermionic degrees of freedom in specific models H. Some 
further parastatistical generalization of the latter interpretation was suggested in Ref. 0. Using 
this method, we shall describe a set of self-similar one-dimensional potentials whose discrete 
spectra are composed from a number of geometric series. Among the corresponding nonlinear 
spectrum generating algebras we shall find a g-deformed Heisenberg-Weyl algebra and a quantum 
conformal algebra su q (l, 1), where parameter q will have the meaning of interaction constant. 

Quantization of the spectral parameter A emerges due to particular boundary conditions im- 
posed upon the wave functions ip(x). Just for an illustration we exhibit one of the possible 
(self-adjoint) conditions: 

V>'(xi) = Mfri), *P'(x 2 ) = M>(x 2 ), (2) 

where X\ and x 2 are two different points on the line, and hi, h 2 are arbitrary real constants. Here 
we would like to mention the powerful restriction on asymptotic growth of eigenvalues {A n } for the 
potentials without singularities in [zi,ar 2 ]. Namely, {A n } can not grow faster than n 2 at n — > oo 
H . This observation immediately restricts the region of applications of g-deformed commutation 
relations which often formally produce exponentially growing discrete spectra. Below we imply 
the condition of square normalizability of wave functions, ip(x) G L 2 [— oo, oo], for the discrete 
energy states. For singular potentials one still will need additional boundary conditions. The 
(quasi)periodic potentials with zonal structure of spectra will also appear in the consideration. 

Let us first present some basic principles of the factorization method. The main object ap- 
pearing after successive factorizations of the stationary one-dimensional Schrodinger operator ([]]) 
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is the following chain of coupled Riccati equations (the dressing chain): 

/;(*) + + - fj+iW = kj = A i+1 - Aj, j e Z (3) 

where kj (Aj) are some constants. The Hamiltonians Lj associated to have the form, 

L 3 =P 2 + ff(x) - f'j{x) + X, = a+aj + Xj, (4) 

af=p±ifj(x), p=—id/dx. 

Conditional isospectrality of Lj and Lj + \ follows from the intertwining relations Lja- = a^Lj+i, 
ajLj = Lj + ia~ . Any exactly solvable spectral problem with infinite number of discrete levels can 
be represented in the form (H), (||) with Aj being the Hamiltonian eigenvalues: 

Ltf® = E$1>V> , E$ — X n+j , n = 0,l,2,... 

^ oc afy£»\ & oc aj^\ (5) 
ip^ix) = exp (- J fj(y)dy). 

In this case A J+ i > Aj and all wave functions ip^ are square integrable, with n nodes inside taken 
coordinate region. For Hamiltonians with finite number of bound states normalizability of "0o^ 
truncates at some j. A large list of potentials, whose spectra are easily found by the factorization 
method, is given in the first paper of Ref.0. In general, however, parameters Aj do not coincide 
with physical eigenvalues since for a given potential they could be chosen arbitrarily. 

In order to solve underdetermined system (||) one has to impose some closure conditions. At 
this stage it is an art of a researcher to find such an ansatz, which allows to generate infinite 
number of fj and kj by simple recurrence relations from fewer entries. Most of the old known 
examples are generated by the choice fj(x) = a(x)j + b(x) + c(x)/j, where a, b, c are some functions 
determined from setting powers of j in (|3]) equal to zero. Three years ago Shabat and Yamilov 
considered the problem from another point of view [0. They suggested to treat (^J) as some 
infinite-dimensional dynamical system and to find potentials corresponding to finite-dimensional 
truncations of the chain. In particular, it was shown that very simple periodic closure conditions: 

fj+N(x) = fj(x), X j+N = X j7 (6) 

for iV odd lead to the potentials with zonal structure of spectra which are described by (hy- 
per)elliptic functions and are called the finite-gap potentials (parameters Aj represent physical 
eigenvalues when they coincide with the boundaries of gaps). First non-trivial example appears 
at iV = 3 and corresponds to the one-gap Lame equation. 
Just slight modification of (H): 

fj +N (x) = fj(x), k j+N = kj, (7) 

describes essentially more complicated potentials, although one has only one additional free pa- 
rameter uj = kj + kj + i + . . . + kj + N^i = Xj + N — Aj 7^ 0. It is easy to see that (0) at N = 1 gives 



2 



the harmonic oscillator problem, whereas (^) is degenerate. The N = 2 system coincides with the 
general conformal quantum mechanical model ||, 

, ( , I, h - k 2 1 h + k 2 . , . 

Mx) = 2 {± k^h- X + —^ X) - (8) 

Already N = 3 case leads to transcendental potentials ||, namely, fj(x) depend now on solutions 
of the Painleve-IV equation |3J : 

= \u>x + f(x), f 2>3 (x) = --/ T^j(f + h), (9) 

f"=Jj + \f + + (f^ + ^ " h)f - M (10) 

To the author's knowledge this is the first example, when Painleve transcendent appears in a 
quantum mechanical context. 

Let us take two unconstrained Hamiltonians Lj and Lj + ^ from the chain (|4]) and assume 
that both are self-adjoint. The map from unnormalized wave functions of Lj+jv to those of Lj is 
performed by successive action of the operators Oj , 

1$(x) = K+^ +N \x), Kf = afaf +1 . . . af +N _ l7 Kj = (K^ 

I\; /-,. v = LjKf, L 3+N KJ = KjLj. (11) 

For a discrete spectrum the labelling of levels differs by some integer, n! = n + m, \m\ < N. Note 
that for superpotentials fj with strong enough singularities some of the Hamiltonians Lj are not 
isospectral because the operators af are not well defined for a subset of physical states. 

Now we restrict potential of Lj+jv to be equal to that of Lj up to some simple transformation. 
Then map eigenstates of Lj onto each other, ip^/ <-> tp^ (there may be zero modes). The 
simplest case is realized when all states are mapped onto themselves, i.e. when n' = n. This 
is equivalent to Lj+jv = Lj, or @. From (|TTD , the operators Kf are seen to be integrals of 
motion |7J. For n'/na spectrum generating situation arises. Taking Xj+N 7^ Aj and substituting 
Lj +N = Lj + Xj+N ~ Xj into ([□]) one gets the standard ladder relations. Integrability of these 
periodic truncations of the chain (H) is analyzed in Refs.[|7|, ^j. Any further generalizations of the 
results mentioned above would involve solvable potentials with infinite number of levels (gaps). 
Let us show that the notion of g-deformation naturally emerges on this track. 



A peculiar potential was found in Ref.[|10f by the following self-similarity constraint: 

f j (x) = qif(qix), k 3 = q 2 ^k, < q < 1, k > 0, (12) 
which gives a solution of @ provided f(x) satisfies the equation 

f{x) + qf(qx) + fix) - q 2 fiqx) = k. (13) 
Quantum algebraic content of this model was uncovered in Ref . 1 1|| , where it was shown that 



g-deformed Heisenberg-Weyl algebra [12 



A~A+ - q 2 A + A- = k, (14) 
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is implemented by the choice 



A+ = (p + if(x))T q , A- = T-\p - if(x)) = (A + )\ 

TJ{x) = Jqf{qx). (15) 

Intertwining relations between A^ and Hamiltonian H = A + A~ — k/(l — q 2 ) easily generate the 
spectrum 

HA ± = q ±2 A ± H, => E n = -kq 2n /{l-q 2 ). (16) 
A deformation of supersymmetric quantum mechanics, inspired by this model, was suggested 



in Ref.[|13|]. The main idea is very simple - one has to replace superpartner Hamiltonian by that 
obtained after affine transformation (i.e. dilatation and translation) and adjust kinetic term to 
the standard form. Degeneracies of levels are removed and energy split is proportional to 1 — q 2 , 
where q is the scaling parameter. Within this scheme, Eq.(|l3l) is a condition of homegeneity of 
magnetic field alond third axis for a spin- 1/2 particle moving on the line. This construction is 
easily generalized to the particular parasupersymmetric model defined by unification of sequential 
members of the chain (^) into diagonal (iV+ 1) x (iV+ 1) matrix. Acting on each subhamiltonian 
by different affine transformation group elements and rearranging kinetic terms one would get 
multiparameter deformation of parasupersymmetric algebraic relations. Following the considera- 
tion of Ref. 0, one may impose various physical restrictions on the matrix Hamiltonian and look 
for explicit form of potentials accepting those constraints. Analyzing such possibilities the author 
have found the following g-periodic closure of the chain (|3|), 

f j+N (x) = qfj(qx), k j+N = q 2 kj. (17) 

It leads to a set of mixed finite-difference-differential equations which describes g-deformation of 
the finite-gap and related potentials discussed in Refs.|7|, |9|. The self-similar system of Ref.[|T0| is 
generated at N — 1. Because of the highly transcendental character of self-similarity and connec- 
tions with the Painleve equations all corresponding potentials may be called as g-transcendental 
ones. 

Let us find a symmetry algebra behind fll?]) . First we rewrite g-periodicity at the Hamiltonian 
level: 

L j+N = q 2 T q L{T~ x + aj, aj = \ j+N - g%, (18) 
where we normalize <jj > 0. Substituting this into ([11]) we get 

LjAf - g 2 ,l; /-, = o r \:. AjL, - g%.l ; = a 5 Aj, (19) 

Aj = Kj~T q , Aj = (Aj) . 

Obtained formulae represent a first part of the quantum algebraic relations determining structure 
of the system. Second part is fixed by a particular g-commutator of A^ following from the 
identities: 

N-l N-l 

A+Aj = J] (Lj - AjAf = II (<?%• + *i - A i+i ). (20) 

i=0 i=0 

As an example, we write one possible equality: 

A+Aj - q 2N A~A+ = VN-tiLj), (21) 
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where Vn-i is a polynomial of the degree N — 1. Formulae (|TP|), ( 2T|) define a particular class 
of nonlinear algebras which may be interpreted as a g-deformation of the polynomial algebra 
of ordinary differential operators discussed in Ref.|| (the latter in turn may be considered as 
the simplest VT-algebras) . Corresponding ladder relations determine the spectrum provided Af 
respect boundary conditions of a problem. The peculiarities of representations of general nonlinear 
algebra [Jo, J±] = ±J±, [J+, J-] = 9 (Jo), where J±,o are some formal operators, were discussed in 
Ref.@. 



For simplicity, we restrict our consideration to the whole line spectral problems with non- 
singular potentials, in which case all a^'s are well defined. Then the index j can always be chosen 
in such a way that Aj will not have zero modes (this will be assumed below). Normalizability of 
physical states is not spoiled by A^-operators, which thus raise and lower energy. As a result, the 
equation Ajip(x) = determines lowest energy state. Suppose that all N independent solutions 
of this equation are normalizable, which corresponds to the ordering Xj < Xj +1 < . . . < A J+ j V _ 1 

and normalizable ipo^ s i n ©■ Energies of all bound states are easily found: 



E U) = a i + J 
1 - g 2 % 



' (Aj - T^)(f r , fo r n = Nr 

(A j+ i - 7^)q 2r , hin = Nr + l 



(22) 



[ fri+N-i ~ 1^2 )g 2r , for n = Nr + N - 1 



By definition < E^ l7 and since q < 1, one has Ej>£ = <Jj/(l — q 2 ), foo( x ) = 0, i.e. the 
potentials are reflectionless, with spectra comprising geometric series. For the continuous 
spectrum the roles of Af are interchanged, action of Aj creates a series of states with exponentially 
growing eigenvalues. Note that at q — >• there will remain only first iV levels corresponding to 
zero value of r in (p2|). The smoothness condition allows to put qfj(qx) = directly in the set of 
equations defining fj(x) and obtain iV-level dressing problem for zero potential. It is known to 
lead to a potential describing fixed time A^-soliton solution of the KdV equation. 

The formula (^) does not work at q > 1. In fact, this region of q may be reached by analytical 
continuation of q < 1 solutions and permutations: 

fi(x, q~ l \ h, k 2 , ■ ■ ■ , k N ) = iqfi(qx/i, q; k N , k N -i, . . . , h), 

f 2 (x,q~ 1 ) =if N (x/i : q), f 3 (x,q~ 1 ) = if N _ x {x/i,q), f N (x, q~ l ) = if 2 (x/i, q), 

where fcj-dependence in the last relations is the same as in the first formula. In Ref.JTT]] the 
N = 1 case was analyzed at q > 1 and the presence of singularities was demonstrated. Analogous 
situation takes place in general - analytic continuation of smooth potentials to imaginary axis 
creates singularities. The latter are moving under the scaling transformation and this breaks 
needed isospectrality of L/s. Note that we do not g-deform Heisenberg equations of motion, so 
that quantum mechanical time evolution coincides with the standard one. However, the time 
evolution of the potentials U(x) as infinite-number soliton solutions of the KdV equation requires 
special consideration. 

In order to elucidate the construction, we consider N = 2 g-periodicity in more detail. Basic 
equations 

f[(x) + &(x) + fl(x) - fl(x) = k u 

ffc) + qfiM + IK*) - q 2 fiM = h, (23) 
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and operators A + — (p + ifi)(p + if2)T q and H = a{a 1 generate the algebra 

HA + — q 2 A + H = (ki + k 2 )A + , A~H — q 2 HA~ = {k x + k 2 )A~, (24) 

A-A + - q A A + A- = q 2 (h(l + q 2 ) + 2k 2 )H + k 2 (h + k 2 ). (25) 

By adding to H some constant one may rewrite equations (|24]) in the form fll6l). Relations 



(|24]), (p5|) at q = 1 define conformal algebra sw(l, 1) and at q ^ 1 they represent a particular 
"quantization" of this Lie algebra (see, e.g., Ref . fl5H ) . 

Let us find fi t2 {x) as formal series near x = 0. Consider first the singular solutions. Permitted 
singularity has a pole character: 



" ^ 2i-l 



Hx) = - + Y,hx 2 -\ /*(*) = -- + £ 

x i=l x i=l 

■ 1 r . r . —h-h- ■ ' ■ a 11 h h - r r • 

C^_^_0^_, g 0j0l _ 3 CjCt _ 

4 2 f-f 2z - 1 + 2a ' y 1 1 ^ 2i - 1 - 2a ' 1 ; 

1 , fci fc 2 , 1 , k 2 _ q 2 k x , 

1_ l-g 2 4 + 2a l-2a j ' Cl ~l-g 2 4-2a l + 2a j ' 
where a is an arbitrary parameter. In general, the series diverge at q — > 1, proper choice of a, 
however, provides the truncated solution (||). In the limit q — ► 0, the function qfi(qx) does not 
vanish, qf±(qx) — > a/a;. Substituting this into fl2"3|) one gets two-level dressing problem for the 
potential oc 1/a; 2 . We guess that at q < 1 there exist such k± j2 that series converge for arbitrarily 
large x. The condition a(a + 1) > 3/4 guarantees that normalizable wave functions and their first 
derivatives vanish at zero 0. The spectrum of such system would arise from only one geometric 
series (second is eliminated by boundary conditions). 

For non-singular at zero solutions one has f± i2 = J2t=o h^'^x 1 , where b^' 2 ^ are two arbitrary 
constants. Again, in general series diverge at q — ► 1. Particular choice of initial conditions gives 
the solution which in this limit corresponds to (§) with a coordinate shift. Depending on the 
values of k± and k 2 the limit q — > recovers either the smooth, two-level potential with bound 
state energies at and k±, or its analytically continued partner. Only first of these corresponds 
to ( p2"D at N = 2, each series belonging to independent representation of su q (l,l). Moreover, 
it is this solution that reduces to the g-oscillator one (0)-(0) (with q replaced by q 1 ^ 2 ) after 
restrictions f 2 (x) = g 1 / 2 /i(g 1 / 2 x), k 2 = qk\. At q — > 1 spectral series become equidistant which 
means that potentials start to be unbounded at space infinities. Because of the nice connection 
with ordinary conformal model ||, we suggest to call the N = 2 system as g-deformed conformal 
quantum mechanics. 

For the iV = 3 system of equations: 

fs(x) + Qfiiqx) + ft{x) - q 2 f 2 x {qx) = k 3 , (27) 
one can exclude f2,z(x) and get 

/ 2 , 3 (z) = -\f{x) =F 2jM(fW + k 2), (28) 
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/(*) = l - ± ^ L h{x) + I ff(y)dy - |(x - (29) 

r (*) = + vmfi^) + KM - \n*) - ft*) -h- \h) - (so) 

where x is a constant of integration and lo = k\ + k 2 + k 3 . At q — 1 one has = /(x) +u;(a; — 
xq)/2 and ( |3~0"1) becomes the PIV equation (|l0f) . So, the system ([29]), (j30T) describes a g-deformation 



of the PIV transcendent |16|. In fact, all functions fi(x) satisfy one combersome equation with 



different choices of the parameters fcj. As a result, the relations ( pq) give new solutions of the 
g-PIV system in terms of a known one: f 2 (x; ki, k 2 , ^3) = fi(x;k 2 ,k 3 ,q 2 ki), f3(x;ki,k 2 ,k- c 
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/i(x; &3, g 2 fci, q 2 k 2 ). Existence of these nonlinear maps is a result of the hidden self-similarity []17] . 



The notion of g-periodicity ( PD and corresponding algebraic relations ([19]), (p0|) constitute 
main results of this paper. However, above we just outlined properties of the self-similar potentials. 
It is quite interesting to know what kind of potentials one gets as a result of deformation of the 
finite-gap potentials, i.e. when for some j one has kj + kj + i + . . . + kj +N _i = 0. There are other 
possibilities in addition to the mentioned ones. For example, the coordinate x and parameter q 
were taken to be real and nothing prevents from the consideration of complex values as well. An 
interesting situation is described when q is a root of unity, q m = 1. From the relations ( |17D one 
easily sees that now 

fj+mN(x) = q m fj(q m x) = fj{x), k j+mN = q 2m kj = kj, (31) 

which is a subcase of @ because 

X j+mN - Xj = (1 - q 2m )(k 3 + k j+1 + ... + k j+N -i)/(l - q 2 ) = 0. 

At N = 1, g 3 = 1 the simplest Lame equation for the equianharmonic Weierstrass function is 
arising [18| . Corresponding spectral problem is known to be solvable [l9j . A charming property of 
the self-similar systems in these cases is that operationally they are naturally characterized not by 
the generators of the order mN polynomial algebras but rather by their m-th operator roots which 
are well defined and satisfy simpler (although unusual) commutation relations. Unfortunately, the 
general analytical structure of the g-transcendents is not known and this does not allow to analyze 
analytical continuation of the solutions found for the specific values of g. Also, there should exist 
some infinite-gap potentials which are reduced to the self-similar ones in the limit of zero widths 
of the gaps. 

Different problems appear when the scaling operator T q (|15|) in the definition of is replaced 
by the translation operator, T a f(x) = f(x + a). Instead of (|13D we then have 

f'{x) + f(x + a) + f{x) - f 2 (x + a) = k. (32) 

This equation may be considered also as a special limiting case of (0). Since q — 1, at k ^ we 
have a realization of the ordinary Heisenberg-Weyl algebra. The full effect of a non-zero parameter 
a in (^) is not known to the author, but purely monotonic solutions are seen to be forbidden. 
The particular solution at k = is given by 

l ^x-xo)-^) 
2 p(x - xq) - p{a) 
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where p(x) is an arbitrary doubly periodic Weierstrass function (note that the values of corre- 
sponding periods do not enter the equation (|32|), which exhibits a non- uniqueness of the solutions). 

Higher dimensional generalizations of the presented construction are not known to the author 
(except of the simple cases when variables separate and the problem becomes effectively one- 
dimensional). Several concluding remarks are in order. First, the factorization method allows 
to replace superpotentials fj(x) by hermitian matrix functions [p0|] , in which case right hand 
sides of Eqs. (||D, (0) are proportional to unit matrices. Second, there may exist an interesting 



interrelation between the described self-similar potentials and the wavelet analysis |21| where 
affine transformations generate orthonormal bases of the Hilbert space. Finally, the splitting of the 
spectrum consisting of one geometric series ([16]) into the one with N terms (|22| ) (at q — 1 there is a 
splitting of the spectrum of a harmonic oscillator, i.e. of the arithmetic progression [[|) illustrates 
the possibility for searching solutions of the dressing chain (|3|) not by expansion in j (or some 
function of j), as it was done in Ref. 0, but by representation j = IN + i — 1, i = 1, . . . , N, I 6 Z 
and expansion in /. 

The results of this paper were partially presented at the XlXth International Colloquium on 
Group Theoretical Methods in Physics (Salamanca, June 1992). The author is indebted to J.- 
M.Lina, A.Shabat, and L.Vinet for valuable discussions and helpful comments. This research was 
supported by the NSERC of Canada. 
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